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ABSTRACT 

In this paper, we present photometry for young star clusters in M31, which are selected from Caldwell 
et al. These star clusters have been observed as part of the Beijing-Arizona-Taiwan-Connecticut 
(BATC) Multicolor Sky Survey from 1995 February to 2008 March. The BATC images including 
these star clusters are taken with 15 intermediate-band filters covering 3000-10000 A. Combined with 
photometry in the GALEX far- and near-ultraviolet, broad-band UBVRI, SDSS ugriz, and infrared 
JHKs of Two Micron All Sky Survey, we obtain their accurate spectral energy distributions (SEDs) 
from 1538 — 20000 A. We derive these star clusters' ages and masses by comparing their SEDs with 
stellar population synthesis models. Our results are in good agreement with previous determinations. 
The mean value of age and mass of young clusters (< 2 Gyr) is about 385 Myr and 2 x 10''^ Mq, 
respectively. There are two distinct peaks in the age distribution, a highest peak at age ^ 60 Myr and 
a secondary peak around 250 Myr, while the mass distribution shows a single peak around 10^ Mq. 
A few young star clusters have two-body relaxation times greater than their ages, indicating that 
those clusters have not been well dynamically relaxed and therefore have not established the thermal 
equilibrium. There are several regions showing aggregations of young star clusters around the 10 kpc 
ring and the outer ring, indicating that the distribution of the young star clusters is well correlated 
with M31's star-forming regions. The young massive star clusters (age < 100 Myr and mass > 10'' Mq) 
show apparent concentration around the ring splitting region, suggesting a recent passage of a satellite 
galaxy (M32) through M31 disk. 

Subject headings: galaxies: individual (M31) - galaxies: young star clusters - galaxies: stellar content 



1. INTRODUCTION 

Star clusters are considered as important tracers for 
understa nding the forma tion and evolution of their host 
galaxies (i San et aril2010[) . Star cluster systems have been 
traditionally separated into two populations-globular 
clusters and open clusters (GCs and OCs)-on their ages, 
masses, metallicities, and positions. However, more re- 
cent studies have discovered that the distinction be- 
tween GCs and O Cs becomes increasingly blurred (see 
iPerina et alll2010l fo r detail s). 

iGasc oigne fc Kronl (|1952f l listed photometric colors 
and magnitudes for star clusters in Magellanic Clouds 
(MCs) and the Fornax dwarf system and divided them 
into two groups. They found that star clusters in blue 
group have central condensation properties similar to 
those of the red group, which were considered as GCs, 
howev er they c ould not be identified with the Galactic 
OCs. 'H odgd (|1961l) termed 23 clusters in the Large 
Magellanic Cloud (LMC)-differing from GCs in their 
relative youth and OCs in their richness and shape-as 
"young populous clusters" , which were called "young 
massive clusters" (YM Cs) or "blue luminou s compact 
clusters" (BLCCs) by Fusi Pecci et al.l (|2005f ). Actually, 
the blue integrated colors for a cluster may be influenced 
by several factors , such as poor metallicity (the lumi- 
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nosity of the horizontal branch), young age (the posi- 
tion of the main-sequence turnoff stars), and some ex- 
otic stellar populations (e.g., blue str agglers, W olf-Rayet 
stars ) . However, several s tudies (e.g.. IWilliams fc Hodgi 
120011 : iBeaslev et al.l I2004D have reached similar conclu- 
sions that the exceedingly blue colors of BLCCs are a di- 
rect c onsequence of their young ages (see lFusi Pecci et all 
[20051 for details). 

M31 is the largest galaxy in the Local Group, and has 
a large number of star clusters, including young clus- 
ters h a ving b een studied by many authors. iBohlin et al.l 
(|1988l I1993D listed 11 objects in M31 classified as 
blue clusters using the UV color s, most of which have 
been proved to be young clusters (TFusi Pecci et al.ll2005l : 
iCaldwell et all[2009l : IPerina et al.ll20'09Ll2010[ ). except for 
B13 3 and B145, which we re stated as a star and an old 
GC (ICaldwell et al.l [20091) . respectively. ICaldwell etahl 
(|2009l 1201 IfT derived ages and masses for a large sam- 
ple of young clusters, and found that these star clus- 
ters are less than 2 Gyr old, and most of them have 
ages between 10® and 10^ yr and masses ranging from 
2.5 X 10^ Mq to 1.5 X 10^ Mq. These authors also 
stated that the young star clusters in M31 show a range 
of structures, most of which have low concentrations typ- 
ical of OCs in the Milky Way (MW), however, there are 
a few with high concent rations similar to the MW GCs. 
iVansevicius et al.l ([20091) carried out a survey of compact 
star clusters in the southwest part of M31, and suggested 
a rich intermediate-mass star cluster population in M31, 
with a typical age range of 30 Myr — 3 Gyr, peaking at 
~ 70 M yr. In order to ascertain the properties of the 
BLCCs. IPerina et all ()2009l . I2010t ) performed an image 
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survey for 20 BLCCs lying in the disk of M31 using the 
Wide Field and Planetary Camera-2 (WFPC2) on the 
Hubble Space Telescope (HST). In addition, another key 
aim of this HST survey was to determine t he fraction of 
contamination of BLCCs by asterisms, since lCohen et ahl 
(I2OO5) suggested that a large fract ion of the putative 
BLCCs may in fact be just asterisms. iCohen et al.l ()2005f ) 
presented the resulting K' images of six very young or 
young star clusters in M31 observed with the Keck laser 
guide star adaptive optics system, and indicated that 
the four youngest out of t hese s ix objects are asterisms. 
However, I Caldwell et al.l (|2009| ) presented a conclusion 
that these four obj ects are true clusters bas ed on spectra. 
The HST images ()Perina et al.|[200l [20I0I ) showed that 
nineteen of the twenty surveyed candidat es are real star 
cluste rs, and one (NB67) is a bright star. iBarmbv et al.l 
()20 09!t) measured surface brightness profiles for 23 bright 
young star clusters using im ages from the WFPC2, in - 
cluding the sample clusters of iPerina et al.l ()2009l |2010| ). 
and derived the structural properties by fitting the sur- 
face brightness profiles to several structural models. The 
authors stated that the sample young clusters are ex- 
pected to dissolve within a few Gyr and will not survive 
to become old GCs, and that young star clusters in M31 
and MCs follow the same fundamental plane relations as 
old GCs of M31, MCs, the MW and NG C 5128, regard- 
less o f their host galaxy environments. iJohnson et al.l 
(j2012l ) presented a M31 stellar cluster catalog utilizing 
the Panchromatic Hubble Andromeda Treasury survey 
data, which will cover ^ 1/3 of M31 disk with multiple 
filters and allow the identification of thousands of star 
clusters. 

The large population of young star clusters reflect a 
high efficiency of cluster formation, possibly triggered by 
a curre nt interaction event between M31 and its satellite 
galaxy ([Gordon et al.ll2006t iBfock et al.l 12006). suggest- 
ing that young star clusters should b e associated with 
the star-forming (SF) regions of M31. iFan et al.l (|2010[ ) 
found that young clusters (< 2 Gyr) are spatially coinci- 
dent with M31's disk, includin g the 10 kpc ring and the 
outer ring (jGordon et al.ll2006l ). Ahhough these authors 
also found the young star clusters in the halo of M31, 
all of the clusters outside of the optical disk of M31 are 
old, globula r clust ers (see IPerina et al.ir2011l for details). 
iKang et ahl ()2012D stated that most of young star clus- 
ters' kin ematics have the thin, rotating disk component 
(see also lRev et al.ll2007l ). The young star clusters' dis- 
tribution has a distinct peak around 10 — 12 kpc from the 
center in M31 disk, and some young star clusters show 
concentration around the 10 kpc ring splitting regions 
near M32 and most of them have systematically younger 
ages (< 100 Myr). Kang et_.al. (2012) also stated that 
the young star clusters show a spatial distribution sim- 
ilar to OB stars, UV SF regions, and dust, all of which 
are important tracers of disk structures. 

Several criteria were developed for selecting young 
clusters from th e int egrated spectrum and colors. 
iFusi Pecci et al.l (|2005[ ) comprehensively studied the 
properties of 67 very blue and likely YMCs in M31 
selected according to their color [{B — V)o < 0.45] 
and/or the str e ngth of spectral index (H/3 > 3.5 A). 
iPeacock et al.l (|2010[ ) presented a catalog of M31 GCs 
based on images from the SDSS and the Wide Field 



CAMera on the United Kingdom Infrared Telescope and 
selected a populat ion of young clust ers with a definition 
of [{g — r)o < 0.3]. IKang et al.l ()2012l ) published a catalog 
of M31 young clusters (< 1 Gyr) and supported the selec- 
tion criteria \(NUV - r)n < 2.5] and [{FUV - r)o < 3.0] 
(jBohlin et all 119931 : iRev et al.l I2007D . These criterions 
may play important roles in distinguishing young from 
old clusters for those whose ages cannot be derived ac- 
curately. 

The formation and disruption of young star clusters 
represent a latter-day e xamp le of the hierarchical for- 
mation of galaxies (Fall 12004 ) . Motivated by that, we 
decided to describe some basic properties of young star 
clusters in M31, such as positions, distributions of ages 
and masses, correlations of the ages and masses with 
structure parameters, which may provide important in- 
formation about the processes involved in their formation 
and disruption. 

In this paper, we will provide photometry of a set of 
young star clusters in M31 using images obtained with 
the Beijing-Arizona-Taiwan-Connecticut (BATC) Mul- 
ticolor Sky Survey Telescope. By comparing the ob- 
served SEDs with the GALEV simple stellar population 
(SSP) models, we derive their ages and masses. This pa- 
per is organized as follows. In Section 2 we present the 
BATC observations of the sample clusters, the relevant 
data-processing steps, and the GALEX (FUV and NUV), 
optical broad-band, SDSS ugriz and 2MASS NIR data 
that are subsequently used in our analysis. In Section 
3 we derive ages and masses of the sample clusters. A 
discussion on the sample young clusters (< 2 Gyr) will 
be given in Section 4. Finally, we will summarize our 
results in Section 5. 



2. SAMPLE OF STAR CLUSTERS, OBSERVATIONS, AND 
DATA REDUCTION 

2.1. Sample of Star Clusters 

Th e sample of star clusters i n this paper is selected 
from lCaldweU et al.l ()2009Ll2011h . who presented a series 
of studies of M31 young and old clusters, respectively;. 
We se l ected 178 young clusters given in Caldwell et al] 
(j2009l I2011D . and fortunately, all the young clusters 
have been observed with 15 intermediate-band filters of 
the BATC photometric system. However, there are 42 
clusters for which we cannot obtain accurate photomet- 
ric measurements with different reasons as following, a) 
some clusters have one or more nearby very bright ob- 
jects; b) some clusters are very close to other objects; 
c) some clusters are very faint and the signal-to-noise 
ratio (SNR) is low; d) some clusters are superimposed 
onto a bright background; e) some clusters are super- 
imposed onto a strongly variable background. In ad- 
dition. There are several remar kable clusters with "ad- 
hered" ()Vansevicius et al.1l2009() objects in our images, 
such as M088 and its neighbor M089, an d G099 and 
C037-G099X (see also iNarbutis et al][2008l) . In a pre- 
vious paper. [Ma et al.l ( 20111 ) presented the SEDs in 15 
intermediate-band filters of the BATC photometric sys- 
tem for one YMC VDB0-B195D and determined its age 
and mass by comparing its SEDs with the theoretical 
evolutionary population synthesis models. Thus, here 
we will analyze the multicolor photometric properties of 
the remaining 135 clusters. 
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Fig. 1. — Spatial distribution of 135 sample star clusters indicated 
by black dots. A box represents a field view of 58' X 5 8'. The 
larg e ellipse is tiie M31 disk/halo boundary as defined by [Racing 
lfT9 91): the two small ellipses are the _D25 isophotes of NGC 205 
(northwest) and M32 (southeast). 

2.2. Archival Images of the BATC Sky Survey for M31 

Field 

The M31 field is part of a galaxy calibration program of 
BATC Multicolor Sky Survey. The BATC program uses 
the 60/90 cm Schmidt Telescope at the Xinglong Sta- 
tion of the National Astronomical Observatories, Chinese 
Academy of Sciences (NAOC). This system includes 15 
intermediate-band filters, covering a range of wavelength 
from 3000 to 10000 A (see lFan et aIlll996L for details). 
Before February 2006, a Ford Aerospace 2k x 2k thick 
CCD camera was applied, which has a pixel size of 15 
/im and a field of view of 58' x 58', resulting in a resolu- 
tion of 1".67 pixerV After February 2006, a new 4k x 4k 
CCD with a pixel size o f 12 /im was use d, with a reso- 
lution of 1".36 pixel-i (|Fan et al.l 12009') . We obtained 
143.9 hours of imaging of the M31 field covering about 
6 square degrees, consisting of 447 images, through the 
set of 15 filters in five observing runs fro m 1995 to 2008 , 
spanning 13 years (see lFan et al.l l2009: W ang et aIll2010L 
for details). 

Figure 1 shows the spatial distribution of the sample 
clusters and the M31 fields observed with the BATC mul- 
ticolor system, in which a box only indicates a field view 
of 58' X 58' of the thick CCD camera. All the sample 
star clusters are indicated with bla ck dot s, with high- 
accuracy coordinates from Caldwell et al.l (|2b09, 201l!), 
which are based on the i mages from the Local Group 
Galaxies Survey (jMassevI [2006i) and the Digitized Sky 
Survey. 

2.3. Integrated Photometry of the Sample Star Clusters 

We processed all the CCD images to apply stan- 
dard procedures including bias subtraction and flat- 
fielding using an automatic data reduction software 
named PIPELINE I, developed for the BATC Muhicolor 
Sky Survey of the CCD images. BATC magnitudes are 



defined and obtained in a similar wa y as for the spec - 
trophotometric AB magnitude system (|Ma et al.ll2009bl) . 
In order to improve the image quality, multiple images 
of the same filter were combined to one, on which the 
magnitudes of the sample star clusters were determined. 
The absolute fiux of the combined images in the central 
field of M31 (M31-1 in Figure 1) was calibrated using 
observations of standard stars, while the absolute flux of 
the combined images of the M31-2 to M31-7 fields was 
calibrated based on secondary standard transformations 
using the M31-1 field (see Fan et al. 2009, for details). 

We performed standard aperture photometry of our 
sample object s using the PHOT routine in DAOPHOT 
(|Stetsonlll987l) . To ensure that we adopted the most ap- 
propriate photometric radius that included all light from 
the objects, we used 9 different aperture sizes (with radii 
of rap = 2.0, 2.5, 3.0, 3.5, 4.0, 4.5,5.0,5.5,6.0 pixels for the 
old CCD, while for the new CCD, the radii were given 
with pixels corresponding to the same arcsecs on the old 
CCD) to determine the magnitude. We also checked the 
aperture radii carefully on the images by visual exami- 
nation in order not to include the light from extraneous 
objects. The local sky background was measured in an 
annulus with an inner radius of rap -1-1.0 pixel and a 
width of 4.0 pixels for the old CCD, and with an inner 
radius of rap -I- 1.0 pixel and a width of 5.0 pixels for the 
new CCD, respectively. 

There are 40 clusters in this paper which are in com- 
mon with our series of previous p apers (Jiang et al 
200l IMa et al.|[200l I2009bl: iFan et al..„200 9: W ang et al 



2O1O0, and the photometric data of these 40 clusters were 
also derived by these studies. We found that most pho- 
tometric data obtained here are in good agreement with 
those obtained by previous studies. We checked the im- 
ages of those clusters with photometric discrepancy, and 
found most of them were loaded near the bulge (B091) or 
in some disk regions with bright or variable background 
(e.g., B210, M020, M023). The different choice of the 
aperture for photometry and the annulus for background 
caused the discrepancy. 

The SEDs for the sample clusters in M31 are listed 
in Table 1. Columns (1) gives the cluster names. 
Columns (2) to (16) present the magnitudes in the 15 
BATC passbands. The la magnitude uncertainties from 
DAOPHOT are listed for each object on the second line 
for the corresponding passbands. For some objects, the 
magnitudes in some filters could not be obtained be- 
cause of low SNR. We should remind that magnitudes 
with an uncertainty larger than 0.3 will not be used in 
the following analysis, although they are listed in Table 
1. Columns (17) is the photometric aperture adopted in 
this paper. 

2.4. GALEX UV, Optical Broad-band, SDSS, and 
2MASS NIR Photometry 

As our series of papers has pointed out, accurate and 
numerous phot ometric points can derive accurate age s 
of star cluster s (Ide Griis et al.l l2003l: lAnders et aLll2004D . 
iKaviraj et al.l (|2007[ ) stated that the UV photometry is 
powerful for age estimation of young stellar populations, 
and the combination of UV photometry with optical ob- 
servati ons enables one to br eak the age-me t allicit y degen- 
eracy. Ide JongI ()1996D and lAnders et all (|2004D showed 
that the age-metallicity degeneracy can be partially bro- 
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ken by adding N IR photometry to optical colors (see 
iMa et all l2009bl and references th e rein). Several pre- 
vious studies (II3armby et alj 120 00^ 'Gal leti et ail 120041: 
iRev et alJ [20071 : iPeacock et alJ 12010: K ang et all \20m 
have provided magnitudes for star clusters in different 
passbands, which will be used to estimate ages of sample 
star clusters in this paper. 

The latest Revised Bologna Catalo gue of M31 GCs and 
cand i dates (hereafter RBC v.4) (Gal leti et al.ll200ll2006l 
[20071 (2009) includes GALEX (FUV and NUV) fluxes 
from[Rcy ct al. (2007), optical broad-band, 2MASS NIR 
magnitudes for 2045 objects. For UBVRI magnitudes 
given in RBC v.4, the relevant photometric uncertain- 
ties are not listed. There fore, we adopted the original 
UBVRI measurements of i Barmbv et al.l (|2000l) as our 
preferred reference, including their published photomet- 
ric errors. For the remaining objects, the UBVRI mag- 
nitudes from RBC v.4 were ad opted, with th e phot omet- 
ric uncertainties set following iGalleti et al.l (I2004D. i.e., 
±0.08 mag in U and ±0.05 mag in BVRI (|Ma et alJ 
I2009bf ). 

In RBC v.4, the 2MASS JHKs magnitudes were trans- 
formed to GIT photometric system (Galleti ct al. 2004). 
However, we needed the original 2MASS JHKs data 
to compare the observed SEDs with the SSP models, 
so we re versed the t ransfo rmation using the equations 
given by iGarpente'^ (|200l . There were no magnitude 
errors for JHKg bands in RBC v.4, and we obtained 
them by compa ring t he photometric data with Figure 
2 of [Carpenterl ()2001h . in which the photometric error 
was shown as a function of magnitude for star s brighter 
than their observational completeness limits ()Ma et al.l 
I2009bl: iWang et al.ll2010[ ). In addition, since RBC v.4 
provided JHKs magnitudes only for a small number of 
sample star clusters, we adopted JHKs magnitudes from 
the 2MASS-6X-PSG catalog, with 6 tim es the normal 
expo sure of 7.2 s on most fields of M31 (jNantais et al.l 
[20061) . There are 3 kinds of magnitudes given by 2MASS- 
6X-PSC catalog, the "default" magnitude, the r = 4" 
aperture and r = 10" aperture magnitudes. We found 
that the r — 4" aperture magnitudes agree well with the 
magnitudes in RBC v.4, while the other 2 kinds of mag- 
nitudes have large discrepancy with magnitudes in RBC 
v.4. However, for the r = 4" aperture magnitudes, when 
JHKs magnitudes are fainter than m = 16 mag, the dis- 
persions are considerable. In this paper, we preferentially 
adopted 2MASS JHKs magnitudes in RBC v.4. For the 
remaining star clusters, we adopted the r = 4" aper- 
ture magnitudes in the 2MASS-6X-PSC catalog when 
t he magnitudes b righte r than m — 16 mag. 

IPeacock et all ([20lol) performed SDSS ugriz photom- 
etry for 1595 M31 clusters and cluster candidates using 
the program SExtractor on drift scan images of M31 ob- 
tained by the SDSS 2.5-m telescope, which are on the 
AB photometric system (see lBertin fc Arnoutslll996l for 
details). We found that there was a magnitude offset 
(> 0.5 mag) between SDSS and UBVRI for 35 objects 
by transforming between the ugriz and UBVR I hand s 
following the transformations from jester et al.l (|2005l ). 
Because the SDSS ugriz magnitudes provided a more ho- 
mogeneous set of photometric measurements, we adopted 
the SDSS magnitudes and abandoned UBVRI magni- 
tudes for these objects in the following analysis. These 
35 objects are flagged with small "a" in Column 1 of 



T able 2. 

iKang et all (|2012[ ) presented a catalog of 700 con- 
firmed star clusters in M31, providing the most extensive 
and updated UV integrated photometry on the AB pho- 
tometric system based on GALEX ima ging, su persed- 
ing the UV photometry published by Rev et al.l ([2007.) , 
which were included by RBC v.4. Therefo re, we used 
themagnitudes of the FUV and NUV from IKang et al.l 
(|2012D as the UV photometry in our following SED fit- 
ting process. 

We hsted the GALEX, optical broad-band, SDSS 
ugriz, and 2MASS NIR photometry of the sample clus- 
ters in Table 2 (Columns 2 to 16), where the photomet- 
ric errors are listed for each object on the second line for 
the corresponding passbands. As we discussed above, the 
magnitudes with an uncertainty larger than 0.3 will not 
be used in the following analysis. 

2.5. Comparison with Previously Published Photometry 

To check our photometry, we transformed the BATG 
intermediate-band system to the broad-band system us- 
in g the relatio n ships between these two systems derived 



m g the relatio n ships t 
bv lZhou et all (|200l : 



S = TOd + 0.2201(TOc- me) +0.1278 ±0.076 and (1) 



V = mg+ 0.3292(m / - m^) + 0.0476 ± 0.027. (2) 

_B-band photometry can be derived from the BATG 
c, d, and e bands, while T^-band magnitude can be ob- 
tained from the BATG f,g, and h bands. Figure 2 
shows a comparison of the B and V photometry of our 
M31 sample objects with previous measurements frorn . 
IBarmbv et all (|2000l ) (circles) and IGalleti et all (pOOl 
(triangles). 

There are several objects with larger offsets (Am > 
0.5), shown with black solid marks in Figure 2 (M045 
in the top panel; B200D, SK036A, and SK068A in the 
bottom panel). The SNRs of M045, B200D, and SK036A 
are low, and SK068A is superimposed onto a bright back- 
ground, thus we cannot derive accurate photometries for 
these four star clusters. 

The mean B and V magnitude differences-in the sense 
of this paper minus others-are (AS) = 0.002±0.213 mag 
and (AV^) = 0.081 ± 0.209, i.e., there is no system offset 
between our magnitudes and previous determinations. 

2.6. Reddening Values 

We required independently determined reddening val- 
ues to estimate ages of the s ample clusters rob ustly 
and accurately . Here we used IKang et all (I2012D and 
IGaldwell et all ([2009h as our reference. IKang et al.l 
(I2012D derived reddening values from three ways: 
1) mean reddening values from ava il able literature 
(IBarmbv et al.l [2000: iFan et al.l [20081: IGaldwell et all 
120091 120111 ): 2) median reddening values of star clus- 
ters located within an annulus at each 2 kpc radius 
from the center of M31, for these star clusters there 
were no available reddening values in the literature; 3) 
for star clusters at distances larger than 22 kpc from 
the center of M31, the foreground reddening value of 
E{B — V)= 0.13 was adopted. Because all of our sample 
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Fig. 2. — Comparison of our newly obtained star cluster photom- 
etry with previous measurements bv lBarmbv et aP Ij2000l ) (circles) 
and Gallcti ct al. (2004) (triangles). The dashed lines enclose ±0.3 
mag in B and V. The black filled circles and tri angles indicate 
the ob jects with photornetry o ffset > 0.5 mag with lBarmbv et al.l 
1I2OOOI) and lGalleti et al.l | |2CI0D . respectively. 

clusters have a projected galactocentric radius smaller 
than 22 kpc, the reddening values for them were de- 
rived from the first two methods. The redde ning value of 
cluste r LGS04131.1_404612 was not given byl Kang et"all 
(120121) ■ so we adop t ed th e value E{B - V) ^ 0.20 
from iCaldwell et aD (|2009| ). Its reddening uncertainty 
was simply adopted half of the reddening value, i.e., 
<7e{b-v) = 0.10. We noticed that, for star cluster 
B449. the reddening value determined by iKang et al.l 
(l2012f) is very differ ent with the value determined by 
ICaldwell et all (I20TI {AE{B - y) = 1.14). We treated 
both age and reddening value as free parameters, and 
determined the reddening value t o be E(B — V) = 0.1 
which was in good agreement with lCaldwell et all ()2011l ). 
So, in this paper ,w e adopted E{B — V) = 0.13 in 
ICaldwell et al.l (|2011l ). The reddening uncertainty for 
B449 was adopted to be 0.07. Column 4 of Table 4 
lists the reddening values adopted for the sample clusters, 
while Column 5 lists the methods for deriving the red- 
dening values (flag = 1 and 2 indicate that the reddening 
values were obtained by the first and second method in 
Kang et al. 2012, respectively; flag = 3 indicates that the 
reddening values were from Caldwell et al. 2009, 2011, 
only for LGS04131. 1.404612 and B449). 

3. AGE AND MASS DETERMINATION 

3.1. Stellar Populations and Synthetic Photometry 

In order to determine the ages and masses of the 
sample star clusters, we compared their SEDs with 
theoretical stellar population synthesis (SPS) mod- 
els. The SSP models of GALEV (e.g_., iKurth et al 



1999; Sch ulz et all 2002: Anders & Fritzc -v. Alvensleben 
200a) were adopted (.Ma et al.,, 2009b : Wang et al.ll2010D 



in this paper, which are based on the Padova isochrones 
(with the most recent versions using the updated Bertelli 
1994 isochrones, including the thermally-pulsin g asymp- 
totic giant-branch [TP-AGB] phase), and a iSalpeteil 
(|1955f ) stellar initial mass function (IMF) with a lower- 
mass limit of 0.10 Mq and the upper-mass limit between 
50 and 70 Mq depending on metallicity. The full set of 
models span the wavelength range from 91 A to 160 /Ltm. 



These models cover ages from 4 x 10^ to 1.6 x 10^° yr, with 
an age resolution of 4 Myr for ages up to 2.35 Gyr, and 20 
Myr for greater ages. The GALEV SSP models include five 
initial metalhcities, Z = 0.0004, 0.004, 0.008, 0.02 (solar 
metallicity), and 0.05. 

Since our observational data were integrated luminosi- 
ties through our set of filters, we convolved the GALEV 
SSP SEDs with the GALEX FUV and NUV, broad-band 
UBVRI, SDSS ugriz, BATC, and 2MASS JHIu fil- 
ter response curves to obtain synthetic ultraviolet, opti- 
cal, and NIR photometry for comparison. The synthetic 
ith filter magnitude in the AB magnitude system can be 
computed as 



= -2.5 log 



- 48.60, 



(3) 



where is the theoretical SED and ipi is the response 
curve of the ith filter of the corresponding photometric 
systems. Here, F^ varies with age and metallicity. 

3.2. Fits 

We used a minimization test to determine which 
GALEV SSP models are most compatible with the ob- 
served SEDs, following 



i=l 



(77 



(4) 



i^°'^{t) is the integrated magnitude in the ith 
filter of a theoretical SSP at metallicity Z and age i, n is 
the number of the filters used for fitting, toJ^"''^ represents 
the intrinsic integrated magnitude in the same filter and 



,.2 

obs.i 



(5) 



Here, aohs,i is the observational uncertainty, and (Tmod,! is 
the uncer tainty associated wit h the model itself, for the 
ith filter. Cha riot et al.l ()1996[ ) estimated the uncertainty 
associated with the term (Tmod.i by comparing the colors 
obtained from different stella r evolutio n ary t r acks an d 
spectral lib rari es. Following iMa et al.l (|2007l l2009allbl . 
I20111I2012D and lWang et~aLl (|2010[ ). weadopted a^^^ , = 
0.05 mag in this paper, cr-ccd is the uncertainty in the 
reddening value, and R\, = A\ J E(B — V), where A\. is 
taken from lCardeUi et al.l (|l989[ ). Ry = Av/E{B-V) = 
3.1, and tTmd,i is the uncertainty of the distance modulus, 
which is always 0.07 from (m - M)o = 24.47 ± 0.07 mag 
(jMcConnachie et al."200^ 

Pcrina ct al. (2009, 2010) determined ages for 20 pos- 
sible YMCs in M31 with metallicity as a free parame- 
ter of their fit. Their results showed that most YMCs 
in M31 were best fi tted with solar metallicity model. 
ICaldwell et al.l ()2009D claimed that it seemed likely the 
young star clusters have supersolar abundances. In this 
paper, the GALEV models of solar metallicity {Z — 0.02) 
were used to fit the intrinsic SEDs for all the sample clus- 
ters here. As an example, we presented the fitting for 
some sample clusters in Figure 3. During the fitting, we 
found that, for a small number of clusters, some photo- 
metric data cannot be fitted with any SSP models. We 
therefore did not use these deviating photometric data 
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Fig. 3. — Best-fitting integrated SEDs of the GALEV SSP models 
shown in relation to the intrinsic SEDs for our sample star clusters. 
The photometric data points are represented by the symbols with 
error bars (vertical error bars for uncertainties and horizontal ones 
for the approximate wavelength coverage for each filter). Open 
circles represent the calculated magnitudes of the model SEDs for 
each filter. 

points to obtain the best fits. These deviating photomet- 
ric data points are the a magnitude of M070, b magnitude 
of M040, k magnitude of V133, m magnitude of M082 
and MlOl, n and p magnitudes of B195, J magnitude 
of B118D, and H and Ks magnitudes of M091. We also 
noticed that, for some star clusters (B091, B305, B319, 
B392, B458, B480, B484, and DA069), the photometries 
in JHKs bands show obvi ous offsets from G ALEX FUV 
and NUV bands (see also lKang et al.|[2012[ ). Consider- 
ing that the UV photometry is po werful for age estima - 
tion of young stellar populations (|Kavirai et alJ IToOTf ). 
we would adopt the GALEX FUV and NUV photome- 
tries in the SED fitting, and abandon the JHKg magni- 
tudes for these clusters. 

The masses of the sample star clusters were determined 
sequentially. The GALEV models provide absolute magni- 
tudes in 77 fihers for SSPs of 10^ M^, includi ng 66 fihers 
of t he HSZ Johns on UBVRI (Landoh 1983), Cousins 
RI (|LandoltJ [l983h . and JHK (BesseU fc Brett 1988f) 
systems. The difference between the intrinsic and model 
absolute magnitudes provides a direct i neasurement of 
the cluster mass, in units of 10^ Mq (see lMa et al.|[20Tll 
for details). We transformed the 2MA SS JHK^ mag- 
nitude s to the photometric system of iBessell fc Bretti 
(|T988l) using the equations given by ICarpenteii (j200li r 
and estimated masses of the clusters using magnitudes 
in all of the U BVRI and JHKs bands. The masses of 
clusters obtained based on the magnitudes in different 
filters were different, therefore, we averaged them as the 



final cluster mass. 

The masses of 22 clusters were not derived, because the 
UBVRIJHKs magnitudes cannot be used: (1) there are 
no UBVRIJHKs magnitudes (LGS04131.1_404612); (2) 
there are no JHKg magnitudes and the UBVRI mag- 
nitudes were abandoned because of the discrepancy with 
the SDSS ugriz magnitudes (e.g., B195); (3) the UBVRI 
magnitudes were abandoned because of the discrepancy 
with the SDSS ugriz magnitudes and the JHKg magni- 
tudes were abandoned because of the discrepancy with 
the GALEX VY magnitudes (e.g., B319). The ages and 
masses of the sample clusters obtained in this paper are 
listed in Table 4. 

3.3. Comparison with Previous Determinations 

In this paper, we determined ages and masses for 135 
star clusters by comparing their multicolor photometries 
with theoretical SP S models. T hese star clusters were 
from lCaldwell et all ()20Q9l |2QTiI ). who presented a series 
of studies of M31 young and old clusters, respectively. 

As discussed in Section 2.3, there are 40 clusters in this 
paper which are in common with our series of previous 
papers, and the ages of 27 cluste rs were also derived by 
our series of st udies (Ji ang et ah 2003; Fan et all [20061 : 
Ma et al.|[200li2009b:.Wang et al.,i 2010). In the study of 
Jiang et al.r(|2003D . the SSP models of Bruzual & Char- 
lot (G. Bruzual & Chariot 1996, unpublis hed) we re used; 
in the studies of Fa n et al. (2006) and Ma et all ([2006), 
the SSP models of Bruzual fc CharlotI ()2003l ) were used. 
In addition, in these three studies, we used only the 
BATC photometries in 13 passbands, and did not in- 
clude UV data, which is powerful to ol for age estima- 
tion of young stellar populations (see iKang et all 120121 
and reference therein). So, we re-estimated the ages for 
these young star clusters with more photometric date 
including UV da ta, a nd with the s ame S SP models as 
iMa et al.l (l2009bD a i id iWang et all (IMol) us ed. In the 
studies of IMa et al.l (l2009bD and lWang et all (I2010D. we 
use d the metall icities obtained bv iBarmbv et aLH 20000 
and lPerrett et a l. (2002) when estimating the ages o f star 
clusters. The metallicities in IBarmbv et al.l ()2000f ) and 
iPerrett et afl ([2002 ) were determined from the Lick in- 
dices whi ch were calibrated from the Galactic old GCs. 
However. lFusi Pecci et al.l (|2005D claimed that young star 
clusters are probably not so metal-poor as d e duced from 
the metallicities obtained by 'Perret t et al.l (|2002( ). and 
concluded that G— band line strength t ends to under- 
estimate [Fe/H] val ues in IPerrett et al.l (|2002[ ) by more 
than 1 dex (see also IKang et a.l.lj2012D~ As a result, most 
of ages obtained in IMa et all (|2009b[ ) and iWang et al.l 
([2OIOD are older than those obtained in this paper. So, 
we re-estimated the a ges for the young star clusters in 
IMa et al.l (|2009bfl and lWang eraP (|2010( ) with the solar 
metallicity SSP models. 

In addition, nine clusters (B476, BHll, M026, M040, 
M045, M053, M057, M058, and M070) were estimated to 
be older than 2 Gyr in th is paper, which are co nsidered to 
be old star clusters (e.g.. lCaldwell et al.ll2009| ). Since this 
paper focused on young clusters, we would not consider 
these nine old clusters in the following analysis. However, 
we pointed out that the ages of six of these nine star 
clusters were estimated to be younger than 2 Gyr by 
previous studies: the age of B476 was estimated to be 1.2 
Gyr bv .CaldweU et al.. (,2009. ) (however, 7.1 Gyr by Fan 
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2010); the age of BHll was estimated to be 1.6 Gyr by 
IVansevicius et all (I2009D; the ag e of M0 40 was estimated 
to be ~ 130 Myr by ' Fan et all ()2010[ ) and ~ 320 Myr 
bv lKang et allj 2012'): the ag e of M 053 was estimated to 
be 1 Gyr bylCaldwell ei~all (12009^ and - 140 Myr by 
iFan et all (|201% the age of M0 58 was estimated to be 
~ 160 Myr by lFan"eral] (120101): the age o f M070 was 
estimated to be 1.2 Gyr by iCaldwell et al.l (|2009D and 
~ 810 Myr by lFan et all ([20101). The ages of these nine 
clusters deriyed by preyious papers and here are listed 
in Table 3. There are three clusters with ages of 4 Myr, 
which is the lowest age limit of GALEV models. Clusters 
KHM31-37 an d V133 were estimated slightly older b y 
other authors ([Caldwell et al.l 120091: iKang et all I2012D . 
w hile B196D was estimated slightly younger 2 Myr) 
by lCaldwell etall |2009). 

In Figure 4, the estimated ages for the young (< 2 
Gyr) star clusters in this paper were com pared with 
those from preyious st u dies (e.g., [Bcaslcy e t al.l 120041: 
Vanseyicius et al.l 120091: iFan et al.! 2010; Perin a et ail 



20101: ICaldwell et al.l l2o69. 2011; Kang ct al. 2012). The 



star clusters with ages of 4 Myr obtained in this paper are 
drawn with open squares in Fig ure 4. There are fiye clus- 
ters in common between i3easley et al. (2004 !) and this 
paper. We can see that the ages of iBeasley et al.l (I2004D 
are in good agreement with ours^ iVanseyicius et al.l 
(|2009f) estimated ages of star clusters located in the 
southern disk of M31 with UBVRI SED-fitting. There 
is an obyious offset between their estimated ages and 
ours, which is caused by some large scatters. If the 
three clusters with age s greater than 3 Gyr giyen by 
IVanseyicius et ail ()2009[ ). which are drawn with arrows 
in Figure 4, are not included, th e systemat i c offse t can 
be reduced to be -0.25 Gyr. iFan et al.l (|2010( ) esti- 
mated ages of star clusters in M31 with multi-band 
(UBVRIJHKs) SED-fitting, and their resu lts agree well 
with ours, with a small offset (- 0.02 Gyr). |P erina et aLl 
(|2009l I2010D determined ages of 20 possible YMCs in 
M31 by comparing the obseryed color magnitude dia- 
grams and the isochr o nes o f different metallicities and 
ages of 'Girard i et al.l (|2002t ) , and estimated masses of 
these clusters based on the Marast on's SSP models o f 
solar metallicity and iSalpeteii ()1955[ ) and iKroupal ()2001[ ) 
IMFs and the IR magnitudes in the 2MASS-6X-PSC 
catalog. In general, the ages obtained in this pa- 
per are in good agreement with the determinations by 
I Perina eTlil . (2010), with a small deyiation (^ 0.06 Gyr). 
There is a small systematic offse t (~ —0.10 Gyr) be- 



tween ICaldwell etall (120091 120 111) and this paper i.e. 
the ages obtained by ICaldweh et al.l (120091 boill) are 



larger than the ages obtained here. IPerina et al.l ( 201(il ) 
als o found the sy s tema tic offset as the ages obtained 
by ICaldwell et al.l ("2009) are l arger than t h e age s ob- 
tained by IPerina et al. (2010l ). IPerina et all (|2010t ) sug- 
gested that this offset is caused by the su per-solar metal- 
licity models (Z = 0.04) adopted by ICaldwell et al.l 
(I2009D when they determined the ages of star clusters. 
IKang et al.l (|2012l ) deriyed ages for young clusters by fit- 
ting the multi-band photometry with model SEDs, and 
their results are in good agreement with ours. 

In Figure 5, we compared the masses of clusters ob- 
tained in this paper with those from preyious stud- 



ies (e.g., 
iFan et al.l 



Beasley et a ll 12004); IV anseyicius et al.' '20091; 
2010^ Perina et al.ll201Q: .Caldwell et al..,200a 
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Fig. 4. — Comparison o f the ages ob tained here with t hose ob- 
tained b y previous works: rBeaslej ;_et_al. (2004), Vanscvicius et al 
1 20091) IFan et al.l ll20Tot). IPerina et al. (2010), Caldw ell et al 



1 2559 . [20ril) . and [Kang et al.l II2012I) . In each panel, the mean 



value of the age differences (ours minus other study) is given, with 
its standard deviation (ct). The error bars of ages from each study 
are also shown. 



[20Tlt IKang eTalllMl ). The masses of two (KHM31- 
37 and V133) of the three clusters with ages of 4 Myr, 
which are drawn with open squares in Figure 5, were 
deriyed to be lower than 10'^ Mq. The masses esti- 
mated in thi s paper are iii good a greement w ith those 
estimated by iBeasley etall (I2004D. IPerina et a l. (2010), 
ICaldweh et all (|2009L 120111 ) . and IKang e~l.i (i2012,) 
There is an obyious offset between IVanseyicius et al.) 
(|2009D and this paper, which is mainly caused by some 
scatters. The cluster with the largest discrepancy is 
B335, with a m a ss es timate of 5 x 10^ Mq by 
iVanseyicius et al.l (120091 ) . and 1.3 x lO'^ Mq in this 
paper. When B335 is excluded, the offset can be re- 
duce d to be ( ;^ -0 .6 x 10^ Mq). The masses estimated 
by F an et al.l (|2010f ) are slightly less than those estimated 
here, with an offset of ^ 1.5 x 10* Mq. 

4. DISCUSSION OF YOUNG STAR CLUSTERS 
4.1. Position 

Figure 6 shows the number, ages and masses of young 
star clusters (< 2 Gyr) as a function of projected radius 
from the center of M31, adopted at ao = 00^42'"44^30 
and So = -f 41°16'09".0 (J2000.q ) following lHm;hra et al.l 
(199_1J) and iPerrett et al.l (|2002D . In the top panel, the 
histogram for the radial distribution of young star clus- 
ters shows clearly two peaks at 4 — 7 kpc and 9 — 11 
kpc, while in the middle panel and bottom panel, wide 
age and mass distributions can be seen in these two peak 
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Fig. 5. — Comparison of the masses obtained here with those ob- 
tained b y previous works: [Beaslcv ct al. (2004), Vansevicius et al. 
(|2009i), I Fan et al .| |[ 20Tol'l. iPer ina ct al. (2010), Ca ldwell et al.i 
J2OO9, 2011), and K ang et al.n 2012). In each panel, the mean 
value of the mass differences (ours minus other study) is given, 
with its standard deviation (cr). The error bars of masses from 
each study are also shown. 



r egions. 

iKang et al.l (|2012l ) presented the radial distribution of 
clusters against the distance from the center of M31, and 
found that the young clusters show two peaks around 
10 — 12 kpc and 13 — 14 kpc. They also found that the 
UV SF regions show two distinct peaks: a main peak 
at ^ 16 kpc and a secondary peak around 11 kpc. In 
addition, a small peak at 5 — 8 kpc in the distribution of 
ages of UV SF regions agains t the projected radius (see 
Figure 19 of IKang et ai1l20Tl can be clearly found. We 
argued that the peak at 4 — 7 kpc obtained in this paper 
should be associated with the peak at 5 — 8 kpc for the 
UV SF regions, while the peak at 9 — 11 kpc obtained 
in this paper correla te with the well-known 10 kpc ring 
(|Gordon et al.l[2006l) . 

Figure 7 displays the spatial distribution and radial 
distribution of the M31 young clusters with different age 
bins: (a) t < 0.1 Gyr; (b) 0.1 Gyr < t < 0.4 Gyr; 
(c) 0.4 Gyr < i < 1 Gyr; (d) 1 Gyr < t < 2 Gyr. In 
the top panel, young star clusters in different age ranges 
are drawn with different marks. The inner, solid ellipse 
and the dashed cont our represent the 10 kpc ring and 
the outer ring from IGordon et al.l ([2006) based on in- 
frared observations with the Multiband Imaging Pho- 
tometer for Spitzer (MIPS) instrument on the Spitzer 
Space Telescope, respectively. The 10 kpc ring was drawn 
with a cent er offset from the M31 nucleus by [5'. 5, 3'.0] 
(|Gordon et al...2006, ) with a radius of 44 arcmin (10 kpc). 




Projected Radius [kpc] 

Fig. 6. — {Top panel) Number histogram of young star clusters 
against projected radius. (Middle panel) Age versus projected ra- 
dius for sample young clusters. [Bottom panel) Mass versus pro- 
jected radius for sample young star clusters. The open rectangles 
show the two peaks around 4 — 7 kpc and 9 — 11 kpc of the radial 
distribution for young star clusters. 

There are several regions drawn with open rectangles 
which show aggregations of young star clusters to dif- 
ferent extents. However, we should point out that these 
aggregations of young star clusters may be caused by the 
projection effect becaus e of the inclination of M31 disk. 
iVansevicius et al.l ()2009D noted two clumps of young clus- 
ters, both of which are located in one rectangle (~ —13 
kpc < X < -9 kpc and -3 kpc < F < kpc). The 
star clusters in this study are spatially coincident with 
the disk and the rings, indicating that the distribution 
of the young star clusters correlates with the galaxy's 
SF regions, wh i ch is consistent with previous studies 
(|Fan et al.ll2010l: IKang et al.ll2012D . In the bottom panel, 
the number of young star clusters in different age bins as 
a function of projected radial distance from the M31 cen- 
ter was shown. We can see that clusters younger than 
0.1 Gyr show most obvious aggregation around the 10 
k pc ring. 

IGordon et all ([20061 ) ran a number of numerical sim- 
ulations of the M31-M32 and M31-NGC 205 interac- 
tions, and assumed a passage of M32 through the disk 
of M31 occurring 20 Myr ago, resulting in a burst of 
star formati o n tha t propagates outward through the disk. 
iBlock et al.l (|2006[ ) suggested that M32 and M31 had an 
almost head-on collision about 210 Myr ago, and M32 
passed through M3 1 disk agai n about 110 Myr ago (see 
Figure 2 of iBlock et al.l [2006 ^. which induced two off- 
center rings-an inner ring with projected dimensions of 
~ 1.5 kpc and the 10 kpc ring. Both of the simulations 
recurred the 10 kpc ring and the observed split. 

We divided our sample star clusters younger than 300 
Myr into six groups, and showed the spatial distribution 
for each group in Figure 8. We can see that only star clus- 
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Fig. 7. — Spatial distribution {top panel) and radial distribution 
{bottom panel) of M31 young star clusters with different age bins: 
(a) t < 0.1 Gyr; (b) 0.1 Gyr < i < 0.4 Gyr; (c) 0.4 Gyr < t < 
1 Gyr; (d) 1 Gyr < t < 2 Gyr. The inner, solid ellipse and the 
dashed contour repr esent the 10 kpc ring and the outer ring from 
[Gordon et al.l I I2006 I). while the dotted ellipse is the M31 disk/halo 
boundary as defined by Racine ( 1991). The several small rectangles 
show the clumps of young clusters to the extents. 

ters with ages 50 Myr — 100 Myr appear around the 10 
kpc ring and the ring sphtting region (—9.5 kpc < X < 
-7.5 kpc and -2.5 kpc < F < -0.5 kpc) (|Kang et all 
I20121 ). indicating that 1) the 10 kpc ring may begin to 
form about 100 Myr ago; 2) M32 passed through the 
southern part of M31 disk around 100 Myr and in turn 
resuhed in the spHt in the form of a h ole. This appears to 
be consistent with the prediction bv iBlock et al.l ()2006f ) 
of a second passage of M32 about 110 Myr ago. After 
the second passage, star clusters formed around the split 
for a long period, since there are a number of star clus- 
ters around th e split with ages younger than 50 Myr. 
iDavidge et al.l (2012) reported that the star formation 
rate (SFR) of the M31 disk would be elevated greatly 
and quickly after an encounter event, and it would finally 
drop when the interstellar medium is depleted and dis- 
rupted. However, from Figure 8, we cannot fi nd evidence 
of ra.dial trend of star clu ster ages (see also iKan g et al.l 
[Ma ICaldwell et al.l[2009l) . 

Figure 9 shows the spatial and radial distribution of 
the M31 young star clusters with different mass bins: (a) 
10^ Mq < M < 10^ Mq; (b) 10^ Mrs, < M < 10^ Mq] 
(c) M > 10* Mq. In the top panel, clusters of these 
three groups are drawn with different marks. The bot- 
tom panel presents the number of young clusters in dif- 
ferent mass bins as a function of projected radial distance 
from M31 center, and it shows that young clusters more 
massive than 10^ Mq are most concentrated nearby the 
10 kpc ring. 

4.2. Age and Mass Distribution 



'- - 4 - 



50 Myr^age<100 Myr 



''^ 



100 Myr^age<150 Myr 



150 Myraoge<200 Myr 



200 Myr^age<250 Myr 



250 Myr^age<300 Myr 



X [kpc] 



X [kpc; 



Fig. 8. — Spatial distribution of six groups of M31 young star 
clusters younger than 300 Myr, divided with same age bin of 50 
Myr. The inner, solid ellipse and the dashed contour represent the 
10 kpc ring and the outer ring from Gordon et al. (2006), w hile the 
dotted ellipse is the M31 disk/halo boundary as defined bv Ra cing 
1)19911 ). The small rectangle represents the ring splitting region in 
the southern part of M31 disk. 
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{2006), wh ile the dotted ellipse is the M31 disk/halo boundary as 
defined by IRacind 1)19911 ') . The several small rectangles show the 
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Fig. 10. — Age and mass distribution of the sample young star 
clusters in this paper. The histograms for age and mass are pre- 
sented with gray colors. 

Figure 10 plots the distribution of estimated ages and 
masses for the young star clusters. A prominent corre- 
lation can be seen that mass increases with age. There 
are two distinct peaks in the age histogram: a highest 
peak at age ~ 60 Myr (log age = 7.8) and a secondary 
peak around 250 Myr (log age = 8.4). The mass distribu- 
tion of the young star clusters show a single peak around 
10^ Mq. The mean values of age and mass of young clus- 
ters are about 385 Myr and 2 x 10** Mq, s Ughtly higher 
than the values presented by IKang et al.l (j2012() , which 
are 300 Myr and 10^ Mq, respectively. Most of our young 
clusters have masses ranging from 10^'^ Mq to 10^ Mq, 
which are more massiv e than OCs in the solar neigh- 
borhood (jPiskunov et a l. 2008), but less massive than 
typic al GCs in the MW (.Mclaughlin fc van der Marell 
[2001 . The lack of young clusters more massive than 
10^ Mq is also note d by iVansevicius et al.l ()2009[ ) and 
iCaldwell et al.l (I2009D. possib l y cau sed by a low-average 
SFR of M31 ()Barmbv et al.l 12001 ) or hidden by dust 
clouds in the d isk due to the inclination angle of M31 
(IVansevicius et al . 2009). 

IPortegies et al.l (|2010;) have listed three phases for the 
evolution of a young star cluster: 1) the first few Myr, 
during which the star formation activity is still proceed- 
ing and the star cluster is rich in gas; 2) a subsequent pe- 
riod after the first supernovae (some 3 Myr after forma- 
tion), in which a young cluster is experiencing a serious 
loss of gas and dust, and stellar mass loss plays an im- 
portant role in the cluster evolution; 3) a later stage that 
stellar dynamical processes dominate the cluster evolu- 
tion. The dividing line between phase 2 and phase 3 may 
be anywhere between 100 Myr and 1 Gyr, and most of 
our young clusters are experiencing the phase 2 or phase 
3.^ 

iChernoff fc Shapirol (|1987[ ) presented that after 5 Gyr, 
both mass and galactic loc ation are impor tant evolution- 
ary parameters for GCs. iSpitzeii (|1958l ) discussed the 
destructive effects of encounters of clusters with giant 



molecular clouds (GMCs), and presented that the dis- 
ruption time for a star cluster varies directly with the 
cluster densit y and is about 200 Myr for a mean density 
of 1 Mq/pc^. iSpitzer fc HarmI (|1958l ) also reported that 
two-body relaxation is effective at destroying low-mass 
clusters and this may account for the scarcity of low- mass 
older clusters. Actually, the two-body relaxation and the 
encounters with GMCs are also imp ortant processes that 
lead to young cluster disrupti on (seelCaldwell et"al] 20091 
and references therein), while IPortegies et al.l (|2010 l) pre- 
sented that mass loss due to stellar evolution is the most 
important process in the young cluster dissolution. It is 
evident that star cluster mass is one key p a ramet er in the 
star cluster evolution. iBoutloukos et al.l (|2003[ ) derived 
an empirical relation between the disruption time and 
the initial mass of star clusters in the solar neighbor- 
hood, Small Magellanic Cloud (SMC), M51, and M33. 
iLamers et al.l (j2005i) determined a disruption time of 1.3 
Gyr for a 10^ Mp^ c luster in the solar neighborhood, 
while ICaldwell etall (|2009( ) reported that most of M31 
young clusters would be destroyed in the next Gyr or so, 
and only some massive and dense ones may survive for a 
longer time. 

Several features are shown in Figure 10: 1) there is an 
obvious gap in the age distribution around 100 Myr. 2) 
there are few clusters older than 400 Myr (log age = 8.6) 
with mass lower than lO'' Mq. Although many low- mass 
clusters can be easy to disr upt, this gap may be caused 
by a selection effect. In fact. lJohnson et all (|2012D found 
that the completeness of M31 ground-based sample drops 
precipitously at mF475w > 18 (MF475W > —6.5), which 
is about 2 x 10** Mq. 3) there are few clusters more 
massive than 10^ Mq, which may be caused by a low- 
average SFR of M31 or the hidden by dust clouds in 
the M31 disk as discussed above. 4) there is a gap of 
clusters with very low masses (~ 10'^ Mq) and younger 
than 30 Myr (log age = 7.5) . These cluster s may be to o 
faint to be sample objects of lCaldwell et al.l (|2009ll201lD . 
indicating that our sampl e is not complete in these age 
and mass ranges (see also ICaldwell et al.ll2009l) . 

Figure 11 shows the age distribution in different mass 
intervals (top panel) and mass function in different 
age intervals (bottom panel). The histograms are de- 
rived using a 0.4-dex bin width with different starting 
values. These distributions contain information about 
the formation and di sruption history of star clusters 
(jFall fc Chandaii l2Q12f ). however, the interpretation of 
the empirical distributions of clusters depe nds strongly 
on ho w incompleteness affects the sample ('Giele s et al.l 
120071) . In the top panel, we can see an obvious gap be- 
fore 40 Myr (log age — 7.6), which is caused by a selec- 
tion effect. The age distribution of the clusters does not 
declines monotonically, with an apparent bend around 
200 Myr (log age = 8.3). We argued that this bend near 
200 Myr may be explained as a burst of cluster forma- 
tion, possibly caused by a current interaction event be- 
tween M31 and its satellite galaxy, such as the collision 
between M31 and M32 about 210 Myr ago suggested by 
iBlock et al.l (j2006[) . The two decline trends starting from 
40 Myr and 200 Myr reflect a rapid disruption of clus- 
ters. Vans evicius et all (|2009f ) noted a peak of the cluster 
age distribution at 70 Myr, and suggested an enhanced 
cluster formation episode at that epoch. In the bottom 
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Fig. 11. — Age distribution of the sample young star clusters with 
differen mass intervals [top panel) and mass function with differen 
age intervals {bottom panel). The histograms are derived using a 
0.4-dex bin width with different starting values. 

panel, the gap in the number of clusters in the low-mass 
regions (log mass < 3.5) is apparenth^ due to a sam- 
ple in completeness, but not physical fVanscvicius et ahl 
[20091 ). The initial mass function for star clusters should 
be slightly steeper (jFall fc Chandad[2012) than what is 
shown here becau se of the short l i fetim es of low-mass 
clusters. Recently, iFall &: Chandaii ()2012l ) compared the 
observed age distributions and mass functions of star 
clusters in the MW, MCs, MSB, M51, and Antennae, 
and found that these distributions of clusters are similar 
in different galaxies. However, due to the incompleteness 
of our cluster sample, partly due to the exclusion of clus- 
ters that cannot derive accurate photometry, we would 
not give any empirical formulas of the distributions for 
age and mass. 

4.3. Correlations with Structure Parameters 

In this section, we will discuss the correlations of ages 
and masses with struc ture pararn eters, which are de- 
rived by King-model (iKind 119661) fits for clusters in 
M31 (|Barmbv et all 120021 120071 [20091) . Because the 
sample clusters are younger than 2 Gyr, the struc- 
ture param eters obtained from the bluer filters are pre- 
ferred (see iBarm bv et al.l 120091 in detail) . There are 
four clusters (B315, B3 19, B368, and B374) which have 
been studied twice by iBarmbv et all (|2002l . |2007| ) and 
iBarmbv et all (I2009D. a nd we would use the new results 
in IBarmbv et all (|2009D . 

Figure 12 shows structure parameters as a function 
of age for young clusters in this paper. Some cor- 
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relations can be seen, the concentration c, defined as 
c = log(rt/ro), decreases with age. The trend is largely 
driven by clusters B342 and B368, both of which have 
large c values (3.98 for B342 and 3.87 for B368). Both 
the scale radius ro and projected core radius Rc increase 
with age. Clusters B342 and B368 have very small rg 
and Rc values and are drawn with arrows in Figure 
12 (The values of ro and Rc are ~ 0.014 pc fo r B342, 
whil e are ^ 0.01 1 pc for B368). lElson et al.l (|l989[ ) 
and lElsonI (|1991[ ) discussed the trend for core radius 
against age, and argued that this trend may represent 
real evolution in the structure of clusters as they grow 
old, partially explained by the effect of mass segregation 
([Mackev fc Gilmori 120031) . or dyna mical effects such a s 
heating by black hole (B H) binaries ([Mackev et al.ll2007D . 
[Wilkinson et al.l ([20031 ) also demonstrated that neither 
large differences in primordial binary fraction nor a tidal 
heating due to differences in the cluster orbits could ac- 
count for the observed trend. The best-fit central surface 
brightness mv o shows a decreasing trend with age, and 
IBarmbv et al.l (2009) argued that this trend may be likely 
due to the fading of stellar population and the increase 
of core radius Rc with age. We also see that the central 
mass density p n decreases wi t h age , although the scat- 
ters are great. IBarmbv et al.l ([2009D presented that the 
central mass density shows very little trend with age for 
both the M31 young clusters and young clusters in the 
MCs. There is no obvious correlation between tr^h, the 
two-body relaxation time at the model-projected half- 
mass radius, and age. The dashed line represents the 
region that tr^h equal to age. It can be seen that most 
clusters (except for DA038 and M091) have ages less 
than tr^h, indicating that these young clusters have not 
been well dynamically relaxed. Because two-body en- 
counters can transfer energy between individual stars and 
then impel the system to establish thermal equilibrium 
([Portegies et al.|[2010D . we argue that these young clus- 
ters have not established thermal equilibrium. 

Figure 13 shows structure parameters as a function of 
mass for the sample young clusters. The concentration 
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Fig. 13. — Structure parameters as a function of mass for the 
sample young star clusters in this paper. 

c increases with mass, although the trend is much weak. 
iFall fc Chandad ()20ia presented c plotted against mass 
for clusters in MCs, and found that there was no corre- 
lation between c and mass. Ac t ually, we found that all 
the clusters in iFall fc Chandaii ()2012D have c less than 
2.5, much smaller than the largest value in our sample 
(~ 4). If we do not include the two clusters B342 and 
B368, the correlation for c with mass nearly disappear. 
Both To and Rc increase with mass, however, the trend is 
largely weaken by cluster B327, which has very small val- 
ues of To and Rc, but larger than those of B342 and B368 
which are drawn with arrows in Figure 13. Both the cen- 
tral surface brightness iiv.o and central mass density po 
decrease weakly with mass, while no obvious correlation 
between tr^h and mass can be seen. 

We checked the surface brightness pro f iles o f B327, 
B342, and B368 displayed in lBarmbv et all ()2009[ ). which 
have very small tq and Rc and very large M v.n and pn, and 
found that these core profiles are cuspy. iBarmby et alJ 
(pOOQ^ concluded that the cores of these clusters did not 
appear to be resolved in the HST/WFPC2 images and 
the structural parameters for these clusters would be un- 
certain if the central cluster luminosity is dominated by 
only a few bright stars. However, if these cuspy core pro- 
files are true integrated properties, which may be better 
fitted by a power-law structure model (e.g.,Eersic 196^, 
the three clusters may have been post core-collapse (see 
iTanvir et aLll2012l in detail). 

4.4. Young Massive Clusters 

YMCs are often related to the violent SF episodes trig- 
gered by galaxy coll isions, mergers, and close encounters 
()de Griis &: Parrnen tier 200 1^. However, ba s ed on a sam- 
ple of 21 nearby spirals, iLarsen &: Richtleil (|1999f) found 
that YMCs can exist in a wide variety of host galaxy en- 
vironments, including quiescent galaxies, and that there 
is no correlation between the morphological type of the 
galaxies and their contents of YMCs. YMCs are dense 
aggregates of young stars, which are also expected to 
be the nurseries for many unusual objects, including 
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Fig. 14. — Spatial distribution for YMCs drawn with different 
sizes of the open circles indicating different mass ranges. The inner, 
solid ellipse and the dashed contou r repr esent the 10 kpc ring and 
the outer ring from Gordon et al^ lj200g ) . wh ile the dotted ellipse 
is the M31 disk/halo boundary as defined bv lRaciri3 11199 II ). The 
small rectangle represents the ring splitting region in the southern 
part of M31 disk, and the two filled black triangles represent M32 
and NGC 205. 



exotic stars, binaries, and BHs (Portegi es et al.|[2010D . 
Many stud ies (.Barm bv et al. 2009; Cald well et alj 120091: 
Vanscviciu s et al 12009: Peacock et al. 20101: iPerina et alJ 
12010: .Portegies et al.li201Ql: iMa et al.ii20lll) that focused 
on M31 YMCs have derived remarkable achievements 
in understanding their stellar populations, structure pa- 
rameters, and dynamical properties. 

There are 13 YMCs in our cluster sample with a 
definition of age < 100 Myr and mass > 10^ Mq 
(jPortegies et al.l[201ClD . Figure 14 shows the spatial dis- 
tribution of the 13 YMCs, while different sizes of the open 
circles indicate YMCs in different mass ranges. The rect- 
angle between the 10 kpc ring and the outer ring repre- 
sents the split in the southern part of the M31 disk, and 
the two black filled triangles represent M32 and NGC 
205. It is not surprising to see that most of the YMCs 
gather around the split, indicating that there has been 
a high-level star forma tion activity, which i s consistent 
with previous studies ([Gordon et al.l 120061 : iKang et all 

5. SUMMARY 

In this paper, we determined the ages and masses for 
a sample of M31 star clusters by comparing the multi- 
color photometry with theoretical SPS models. Multi- 
color photometric data are from the GALEX FUV and 
NUV, broadband UBVRI, SDSS ugriz, 15 intermediate- 
band fihers of BATC, and 2MASS JHKs, which consti- 
tute the SEDs covering 1538 - 20000 A. 

We made a discussion on the spatial distribution, dis- 
tribution of ages and masses, correlations of ages and 
masses with structure parameters for the sample young 
clusters (< 2 Gyr). The mean value of age and mass of 
young clusters is about 385 Myr and 2 x 10^ Mq, respec- 
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tively. There are two distinct peaks in the age distribu- 
tion, a highest peak at age ^ 60 Myr and a secondary 
peak around 250 Myr, while the mass distribution shows 
a single peak around 10"* Mq. There are several regions 
showing aggregations of young clusters around the 10 kpc 
ring and the outer ring, indicating that the distribution 
of the young clusters correlates well with M31's SF re- 
gions. The ages and masses show apparent correlations 
with some structure parameters. We also found the cor- 
relation between core radius Rc and age, which has been 
studied by many authors. A few young clusters have the 
two-body relaxation times t^.ft greater than their ages, 
indicating that they have not been well dynamically re- 
laxed. We argued that these young clusters have not 
established the thermal equilibrium. 
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RATn TMTERMBTITaTT^-TiAMn 



TABLE 1 



Object a b c d e 
(mag) (mag) (mag) (mag) (mag) 



fghijkmnop 7?ap 
(mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (") 



B006D-D036 

BOIOD 

B012D-D039 

B014D 

B015D-D041 

B017D 

B018-G071 

B019D 

B035D 

B040-G102 



20.07 
0.142 

20.89 
0.254 



20.22 
0.216 
19.81 
0.274 



21.08 
0.399 
18.68 
0.097 



19.46 
0.121 
20.03 
0.161 
20.14 
0.261 
19.80 
0.158 
19.51 
0.135 
19.44 
0.154 
18.78 
0.052 
19.81 
0.176 
19.74 
0.192 
17.73 
0.020 



18.99 
0.049 
19.79 
0.045 
19.46 
0.058 
19.64 
0.066 
19.26 
0.051 
18.98 
0.042 
18.32 
0.061 
20.18 
0.139 
19.14 
0.032 
17.64 
0.030 



18.76 
0.056 
19.28 
0.038 
19.19 
0.087 
19.40 
0.099 
18.92 
0.067 
18.51 
0.051 
17.93 
0.083 
19.90 
0.179 
18.91 
0.047 
17.65 
0.056 



18.75 
0.066 
19.23 
0.039 
19.14 
0.093 
19.21 
0.080 
18.84 
0.057 
18.32 
0.046 
17.88 
0.048 
19.48 
0.140 
18.67 
0.040 
17.54 
0.036 



18.61 
0.066 
18.99 
0.034 
18.83 
0.080 
19.04 
0.074 
18.61 
0.053 
18.09 
0.042 
17.80 
0.037 
19.23 
0.128 
18.56 
0.040 
17.48 
0.027 



18.51 
0.059 
18.96 
0.048 
18.85 
0.113 
18.79 
0.067 
18.46 
0.059 
17.97 
0.055 
17.45 
0.037 
18.96 
0.124 
18.41 
0.048 
17.37 
0.037 



18.44 
0.063 
18.78 
0.040 
18.73 
0.114 
18.73 
0.070 
18.34 
0.058 
17.87 
0.048 
17.39 
0.033 
18.79 
0.111 
18.35 
0.046 
17.23 
0.027 



18.32 
0.074 
18.68 
0.045 
18.55 
0.111 
18.59 
0.072 
18.03 
0.052 
17.68 
0.057 
17.20 
0.026 
18.51 
0.108 
18.13 
0.053 
17.18 
0.019 



18.63 
0.052 
18.47 
0.119 
18.43 
0.074 
17.98 
0.061 
17.53 
0.059 
17.13 
0.041 
18.44 
0.101 
18.10 
0.054 
17.18 
0.038 



18.38 
0.108 
18.54 
0.072 
18.36 
0.110 
18.33 
0.093 
18.04 
0.077 
17.41 
0.075 
16.84 
0.045 
18.18 
0.099 
18.05 
0.081 
16.94 
0.055 



18.29 
0.115 
18.57 
0.059 
18.34 
0.117 
18.16 
0.067 
17.86 
0.060 
17.22 
0.062 
16.95 
0.034 
18.39 
0.103 
18.01 
0.060 
17.16 
0.040 



18.45 
0.100 
18.25 
0.139 
18.41 
0.128 
17.91 
0.098 



16.79 
0.048 
18.15 
0.130 
18.18 
0.117 
16.98 
0.054 



18.10 
0.130 
18.37 
0.078 
17.98 
0.105 
18.02 
0.097 
17.65 
0.068 
17.14 
0.099 
16.71 
0.081 
17.84 
0.103 
17.94 
0.120 
16.93 
0.112 



18.11 
0.124 
18.01 
0.167 
17.89 
0.148 
17.65 
0.107 
17.32 
0.157 
16.86 
0.125 
17.95 
0.163 
18.08 
0.177 
16.90 
0.124 



5.0 



4.2 



5.0 



5.0 



5.0 



6.7 



5.8 



5.0 



5.0 



5.8 



TABLE 2 

GALEX, OPTICAL BROAD-BAND, SDSS, AND 2MASS NIR PHOTOMETRY OF 135 SAMPLE STAR CLUSTERS IN M31. 



Object 


FUV 


NUV 


U 


B 


V 


R 


I 


u 


9 


r 


i 


z 


J 


H 






(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


B006D-D036 


21.50 


20.87 






18.69 






19.72 


18.44 


18.30 


18.10 


17.99 










0.100 


0.060 






0.050 






0.114 


0.069 


0.077 


0.090 


0.134 








BOIOD 










18.92 






20.80 


19.11 


18.62 


18.40 


18.27 


















0.050 






0.224 


0.114 


0.120 


0.127 


0.171 








B012D-D039 










19.05 






20.60 


19.03 


18.48 


18.13 


17.81 


















0.050 






0.187 


0.052 


0.068 


0.085 


0.123 








B014D 




22.81 






18.53 






20.54 


19.05 


18.49 


18.14 






15.79 


15.72 






0.380 






0.050 






0.131 


0.044 


0.057 


0.072 






0.133 


0.224 


B015D-D041 










18.70 






20.39 


18.70 


18.01 


17.71 


17.42 




15.92 


15.72 












0.050 






0.183 


0.091 


0.098 


0.108 


0.131 




0.265 


0.151 


B017D 










18.23 






20.05 


18.32 


17.77 


17.39 


17.17 


15.21 
















0.050 






0.132 


0.048 


0.057 


0.066 


0.088 


0.094 






B018-G071 




22.19 


18.47 


18.25 


17.53 


17.00 


16.38 


19.39 


17.85 


17.22 


16.82 


16.63 


15.46 


14.78 


14.60 






0.140 


0.080 


0.050 


0.050 


0.050 


0.050 


0.088 


0.059 


0.063 


0.067 


0.076 


0.094 


0.086 


0.112 


B019D 










18.93 






21.14 


19.38 


18.61 


18.11 


17.59 


16.07 


15.73 


14.79 












0.050 






0.246 


0.043 


0.064 


0.082 


0.110 


0.112 


0.112 


0.112 


B035D 


22.37 


22.95 






18.48 






20.05 


18.68 


18.27 


17.91 


17.58 


15.99 


15.35 


15.18 




0.340 


0.310 






0.050 






0.162 


0.111 


0.121 


0.135 


0.164 


0.094 


0.112 


0.130 


B040-G102 


19.84 


19.30 


17.42 


17.69 


17.30 


16.92 


16.63 


18.32 


17.38 


17.19 


17.05 


16.98 


16.07 


14.96 


15.51 




0.030 


0.020 


0.046 


0.022 


0.010 


0.022 


0.022 


0.048 


0.035 


0.039 


0.045 


0.059 


0.112 


0.086 


0.130 
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TABLE 3 

Age comparison for the nine clusters older than 2 Gyr obtained in this paper with previous studies. 



Object 


Log (Age)" 


Log (Age)*" 




(yr) 


(yr) 


B476-D074 


9.08' 9.85^ 


9.79 ± 0.44 


BHll 


9.20" 


9.67 ±0.97 


M026 


9.93^ 


10.11 ±0.23 


M040 


8.11^ 8.50^ 


9.38 ±0.57 


M045 




9.38 ±0.50 


M053 


9.00' 8.16^ 


9.71 ±0.83 


M057 




9.77 ± 1.03 


M058 


8.21^ 


10.09 ±0.29 


M070 


9.08' 8.91^ 


9.97 ±0.65 



"The age estimates are from lCaldwell et al.l l(2009l.[20Tll ') (ref=l'). lFan et al.ll|2010h ('ref=2'). lKang et aLlpOTl l (ref=3), and lVansevicius et al.l 

||2009|) (ref=4). 

''The ages obtained here. 



TABLE 4 

Ages, masses, and reddening values of 135 sample star clusters in M31. 



Object Log (Age) Log (Mass) E{B - V) Flag" 
(yr) (Me) 



B006D-D036 


7.78 


± 


0.06 


3.76 


± 


0.02 


0.33 


± 


0.05 


1 


BOIOD 


8.70 


± 


0.07 


4.13 


± 


0.02 


0.25 


± 


0.09 


1 


B012D-D039 


7.76 


± 


0.08 


3.83 


± 


0.02 


0.52 


± 


0.05 


1 


B014D 


7.75 


± 


0.11 


3.94 


± 


0.07 


0.50 


± 


0.05 


1 


B015D-D041 


7.78 


± 


0.03 


4.02 


± 


0.13 


0.65 


± 


0.05 


1 


B017D 


9.09 


± 


0.07 


4.76 


± 


0.11 


0.23 


± 


0.18 


2 


B018-G071 


8.92 


± 


0.04 


4.76 


± 


0.05 


0.20 


± 


0.06 


1 


B019D 


9.09 


± 


0.02 


4.46 


± 


0.10 


0.30 


± 


0.05 


1 


B035D 


8.60 


± 


0.03 


4.26 


± 


0.05 


0.25 


± 


0.09 


1 


B040-G102 


7.76 


± 


0.05 


4.13 


± 


0.10 


0.25 


± 


0.09 


1 



" The reddening values are from lKang eFal] (goTl) (Flag = 1, 2) and lCaldwell et all Ip009l . [20TI1 ') (Flag = 3). 



